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Lecture  7 
PULSATING  FLAMES 

In  section  5-  it  was  found  that  plane  'IFFs  of  sufficiently  large 
Levis  number  are  unstable.  Since  Im(ot)  ^  C  on  the  stability  boundary,  the 
instability  is  likely  to  result  in  either  a  pulsating  flame  or  a  flame 
that  supports  traveling  waves.  Such  flames  are  the  subject  of  this 
lecture. 

One  difficulty  that  immediately  confronts  us  is  that,  in  contrast 

to  the  ubiquitous  nature  of  cellular  instabilities,  pulsating  instabilities 

are  not  ordinarily  seen.  The  reason  seems  to  be  the  large  values  of  L  are 

2 

needed;  according  to  the  theory^Y^i/T^  must  exceed  32/3  (or  U(l+/3) 
if  the  disturbances  are  one-dimensional).  There  is  evidence  that  fuel- 
rich  hydrogen-bromine  mixtures  night  attain  such  values  since  Golovichev, 
Grishin,  Agranat  and  Bartsun  (1978)  obtained  oscillations  in  a  numerical 
study,  but  there  is  no  similar  evidence  for  more  commonplace  gas  mixtures. 

For  this  reason  we  must  turn  from  the  commonplace  and  deal  either  with 
unusual  combustible  materials  or  else  with  special  configurations,  in  order 
to  uncover  pulsating  flames.  Our  discussion  will  start  with  thermites, 
which  are  solids  that  burn  to  form  solids  (a  phenomenon  that  is 
appropriately  called  gasless  combustion).  There  is  no  significant 
diffusion  of  mass,  so  that  L  is  effectively  infinite. 

1.  Solid  Combustion 

Experiments  by  Merzhanov,  Filonenko,  and  3orovinskaya  on  thermites 
composed  of  niobium  and  boron  revealed  pulsating  instabilities,  as  did 
numerical  solutions  obtained  by  Shkadinsky,  Khaikin,  and  Merzhanov. 

The  latter  examined  the  equations 
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3T/3t  -  32T/3x2  =  -3Y/3t  =  ft  (l) 

where 

ft  =  VYe~9/T  with  V  =  DM2,  (2) 

r 

and  uncovered  a  critical  value  0  of  the  activation  energy:  for  9  <  0 

c  c 

the  propagation  is  steady,  but  for  0  >  0c  only  pulsating  propagation  is 
seen.  The  prediction  of  pulsating  propagation  is  consistent  with  the 
NEF  analysis  in  section  5-  ,  where  oscillatory  instability  was  found  for 

sufficiently  large  L  in  the  linit  0  -*■  »;  but  activation-energy  asymptotics 
has  nothing  to  say  about  a  phenomenon  (here  the  switch  to  steady  propaga¬ 
tion)  occuring  at  some  finite  value  of  0.  Even  though  it  is  not  observed 
either  experimentally  or  numerically  for  large  enough  0 ,  a  steady  wave 
can  nevertheless  be  constructed  by  means  of  activation-energy  asymptotics; 
we  shall  start  out  discussion  by  doing  so. 

The  following  boundary-value  problem  presents  itself  in  a  frame  moving 
with  the  flame  sheet: 

dT/dx  -  d2T/dx2  =  -dY/dx  =  PYe~0/T,  (3) 

T  -►  Tf,  ?  +  T^  as  x  -*•  -® ,  T  bounded ,  Y  +  0  as  x  +  +■.  (h) 

The  solution  outside  the  flame  sheet  is 


T  is  continuous  across  the  flame  sheet,  as  for  the  finite-L  problem  (cf. 
section  2.  ),  but  Y  Jumps.  As  a  consequence,  the  structure  of  the 

reaction  zone  is  somewhat  different. 
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The  structure  variable  is,  as  usual. 


and 


£  =  0x 


(6) 

(7) 


but  we  now  find  the  equations 

o  o  o  -  .  — ®/^j. 

rd'  $/d£  =  -dY  /d£  =  DY  e-9  with  V  =  Vs  /9 


(8) 


2 

for  the  temperature  perturbation  =  -T^d  and  the  leading  term  Yq  in 
the  mass-fraction  expansion-  Since  both  $  and  Y  vanish  as  £  +•, 

we  have 

Yq  =  -T^dd/dS,  (9) 

so  that  only  the  temperature  equation 


d2*/d?2  «  -PCd*/dOe“* 


(10) 


remains.  To  match  with  the  solution  outside  the  flame  sheet,  the  usual 
boundary  conditions 

d  =  -Yf£/?2  +  o ( 1 )  as  £  d  =  o(l)  as  £  +~  (ll) 

must  be  applied.  The  first  integral 

dd/d£  =  CKe'^-l)  (12) 

of  equation  (10)  is  obtained  by  using  the  boundary  condition  (lib);  then 
the  boundary  condition  (lla)  leads  to  the  burning  rate 

-9/2T 

M  =  /D  T.e  /STJ.  (13) 

r  o  j. 
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This  result  has  apparently  not  been  written  down  before,  but  others 
have  derived  it  (e.g.  Peters  1982). 

2.  The  Delta-Function  Model 

We  have  already  noted  that  there  is  little  point  in  investigating 
the  stability  of  this  solution  using  activation-energy  asymptotics. 
Instead,  following  Matkovsky  and  Sivashinsky,  we  shall  use  a  delta- 
function  model  suggested  by  the  asymptotics  above.  Thus,  the  strength 
of  the  delta  function  that  replaces  the  Arrhenius  term  will  be  defined 
so  that  the  mass  flux  through  the  flame  sheet,  in  quite  general  circum¬ 
stances,  is 

-e/2T* 

M  =  /DT*e  /JTjdi  (lU) 

here  is  the  flame  temperature.  The  dimensionless  mass  flux  is  then 

M/Mr  =  (VTb)exp[6(TJt-Tb)/2TbT*],  (15) 

an  expression  that  will  be  simplified  before  use.  For  0  large  (but  not 

necessarily  infinite),  the  preexponential  factor  is  not  significant 

and  can  be  replaced  by  1;  in  addition,  for  small  deviations  of  T*  from 

T.  (such  as  occur  in  a  linear  stability  analysis)  the  exponent  can  be 
0 

p 

replaced  by  0(T--T.  )/2T.  .  The  formula  (15)  then  becomes 

V  =  M/M  =  exp[e(T*-?b)/2?b],  ( 1 6 ) 

since  the  dimensionless  density  may  be  taken  to  be  1;  here  W  is  the 
wave  speed . 

The  Arrhenius  term  (2)  is  now  replaced  by 
n  =  YfWd(n)  with  n  =  x-F(0,0,t), 


(17) 
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vhere  A  is  the  Dirac  delta  function  and 

x  =  ?(y,z,t)  (18) 

is  the  flame  sheet,  written  in  a  (fixed)  coordinate  system  chosen  so  that 
the  x-axis  coincides  with  the  normal  at  the  point  of  interest  at  the  instant 
considered.  The  equations 

3T/3t  -  V2T  =  -3Y/3t  =  n,  (19) 

which  generalize  the  one-dimensional  ones  (l),  then  show  that  the  wave 
speed  -F  (0 ,0  ,t)  is  Just  W  and  that 

w 

6(T)  =0,  5(3T/3n)  =  -YfW,  «(Y)  =  Yf.  (20) 

Precisely  the  same  formulas  can  be  obtained  by  applying  activation- 
energy  asymptotics  to  the  flame  sheet  in  the  unsteady  case  (i.e.  with 
M  i  and  a  flame  temperature  T#  within  0(0  1)  of  T^)  if  it  is 
assumed,  without  Justification,  that  there  is  no  significant  temperature 
gradient  behind  the  flame  sheet.  What  sets  the  delta-function  method 
apart  from  activation-energy  asymptotics  is  that,  in  solving  equations  (19), 
8  is  treated  as  a  finite  parameter,  and  that  there  is  no  requirement  for 
the  outer  solution  to  match  the  inner  in  the  limit  8 

3.  Stability  of  Thermite  Flames. 

Linear  stability  of  the  plane  wave  is  investigated  in  the  manner  of 
section  5.  for  a  plane  XEF.  Disturbances  proportional  to  exp(at+iky) 
are  sought,  resulting  in  the  dispersion  relation 

(2a+e)/l+t(a+k2)  =  G(2a+l)  with  3  =  Yf0/2T2.  (21) 


The  neutral  stability  curve 
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0  =  [ 6k2+2+(2k2+l )  /l6k2+5 ] / ( ^k2+l )  (22) 

is  shown  in  figure  1;  on  it 

a  =  ±i  0s(l+ltk2)2/2  (23) 

is  everywhere  non-zero,  i.e.  the  neutral  modes  are  oscillatory.  Resenblance 
no  the  right  stability  boundary  for  NE7s  in  figure  5.  is  striking.  Note 
that  the  extreme  value 

0  =4.  (2U) 

e 

is  large  enough  to  give  credence  to  the  notion  that  most  of  the  chemical 
reaction  is  confined  to  a  thin  sheet. 

The  results  suggest  that,  for  the  one-dimensional  equations  (l) 
steady  propagation  is  possibly  only  if 

3  <  0  =  2  +  /5;  (25) 

c 

otherwise  pulsating  combustion  occurs.  This  is  in  agreement  with  the 

experimental  and  numerical  results  cited  at  the  beginning  of  section  1. 

Additional  evidence  is  afforded  by  Matkovsky  and  Sivashinsky ' s  demonstration 

that  the  passage  of  0  through  gives  a  supercritical  Hopf  bifurcation. 

Traveling-wave  instabilities,  rather  than  pulsations,  will  occur  if 

disturbances  of  non-zero  wavenumber  are  permitted  by  the  lateral  boundary 

conditions.  The  effect  is  strikingly  seen  for  propagation  through  an 

insulated  circular  cylinder  of  thermite  (Sivashinsky  19Sl).  Now  distur- 

b=rc°s  orooortional  to  elniJ  (kr)  are  sought,  where  r,d  are  oolar 

n 

coordinates  in  the  y,z-plane  and  n  is  a  non-negative  integer.  Once  more 
the  dispersion  relation  (21),  the  neutral  stability  curve  (22),  and  the 
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neutral -mode  frequency  (23)  are  obtained;  but  now  not  an  values  of  k  are 
admissible .  Thermal  insulation  of  the  surface  r  =  a  requires 


J'(ka)  =  0,  i.e.  ka  =  J '  (26) 

n  n,m 

in  a  standard  notation  for  zeros  of  derivatives  of  Bessel  functions.  The 
first  seven  allowable  values  of  ka  are 


J’  =  0,  J*  =  1.8U,  y  .  =  3.05,  j*  =  3.83, 

0,1  1,1  2,1  0,2 


ji  =  U.20,  y  .  =  5.32,  y  =  5-3? 
3*1  ^  l 


(27) 


according  to  Olver 


(196U). 

k 

n,m 


Some  of  the  eigenvalues 


(28) 


are  marked  on  the  neutral  stability  curve  in  figure  1  for  a  =  2  and  k. 

These  two  cases  illustrate  the  general  movement  of  the  k  -points  down 

n,m 

the  curve  as  a  increases. 

Consider  now  what  happens  for  cylinders  of  different  size  as  0 
is  increased  up  to  the  first  onset  of  instability.  (This  is  not  easily 
done  in  a  practical  context.)  Only  the  discrete  points  on  the  neutral 
stability  curve  corresponding  to  the  values  (28)  of  k  are  relevant,  and 
which  mode  will  be  triggered  first  depends  on  the  value  of  a.  As  0 
is  further  increased  the  mode  becomes  an  admissible  instability  that 
developes  a  definite  nonlinear  form  of  the  same  general  character. 

For 

a  <  2.05  (29) 


the  point  1  lies  to  the  right  of  kQ  ±  on  the  neutral  curve;  this  is 
exemplified  by  a  *  2  in  figure  1.  The  first  manifestation  of  instability 
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will  be  plane  pulsations  of  frequency 

w  =±/0c/2;  (30) 

the  corresponding  expressions  for  flame  temperature  and  location  are 

T#  =  T^+ecosut,  x#  =  -t-(8cu/Yf )sin  ut,  (31) 

where  e  is  the  (linear)  disturbance  amplitude.  The  speed  of  the  flame 
is  greater  or  less  than  1  accordingly  as  its  temperature  is  greater  or 
less  than  T.  •  The  temperature  gradient 

sRe[(cos  ut-i  sin  ut )  (l-A+Uiui)  ]/2  (32) 

behind  the  flame  does  not  vanish  but  flucturates  about  zero  (cf .  the  remark 
at  the  end  of  the  last  section). 

The  model  is  based  on  the  assumption  that  the  reactant  is  consumed 
completely,  but  this  is  not  the  case  in  practice.  Indeed,  it  is  sometimes 
possible  to  propagate  a  flame  through  the  same  material  twice.  It  is  to 
be  expected  that  fluctuations  in  temperature  gradient  at  the  reaction  front 
will  result  in  a  layered  burnt  state;  in  the  context  of  activation-energy 
asymptotics,  negative  temperature  gradients  behind  the  reaction  zone  permit 
reactant  leakage.  Merzhanov,  Filonenko,  and  Borovinskaya  noted  layered 
structure  in  burnt  thermites,  with  a  layer  for  each  pulsation. 

A  quite  different  phenomenon  occurs  if  the  first  manifestation  of 
instability  is  associated  with  J '  ,  as  is  the  case  for 

2.05  <  a  <  4.88;  (33) 

r.ov  >  ^  is  the  leftmost  point  on  the  neutral  curve,  as  is  exemplified 

I  ,1 

by  a  =  4  in  figure  1.  The  expressions  (31)  are  replaced  by 
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T„  =  T.o+ecos(ut+!}.)J1(k1  ^r)  ,  x#  =  -t-(8ca>/Yf  Hintwt+ifr).^^  ^r) ,  (3*0 


where  now 


(35) 


The  .  isotherms  of  the  flame  temperature  are  shown  in  figure  2;  these 

.  spin,  in  either  direction,  with  the  frequency  (35),  producing  a  single 
hot-spot  traveling  in  helical  path  on  the  surface  of  the  cylinder.  Such 
hot  spots  were  observed  by  Merzhanov ,  Filonenko  and  Borovinskaya. 

As  the  radius  a  increases  beyond  U.88  for  a  certain  interval, 

£  » ■*« 

becomes  the  leftmost  point  on  the  neutral  curve;  there  are  then  two  spinning 
hot  spots  at  opposite  ends  of  a  diameter.  Hot  spots  can  also  occur  in  the 
interior,  for  example  for  k^  ^ ,  which  is  the  leftmost  point  for  an  interval 
of  still  larger  values  of  a. 

Other  cross  sections  give  rise  to  their  own  distinctive  sets  of  admissible 
wavenumbers  and  isotherm  patterns.  For  rectangles,  judicious  choice  of 
proportion  leads  to  2  or  even  3  modes  simultaneously  characterizing  the 
onset  of  instability.  A  small  change  in  the  proportion  will  cause  the 
corresponding  eigenvalue  to  split,  leading  to  secondary  or  tertiary 
bifurcations.  Matkcwsky  &  Olagunju  (1982)  have  carried  out  the  unimodal 
bifurcation  analysis  for  circular  cross  sections;  Margolis  &  Matkowsky 
(1982)  have  considered  the  multi-modal  analysis  for  rectangular  cross 
sections. 

The  stability  boundary  identified  here,  coupled  with  the  nature  of 
the  instabilities,  shows  that  we  are  dealing  with  the  analog  of  the  right 


stability  boundary  in  the  tFEF  analysis.  This  suggests  that  SVFs,  which 
lie  between  NEFs  (with  L  close  to  l)  and  thermite  flames  (with  £.=•), 
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should  also  exhibit  pulsations  as  their  instability;  however,  for  L  >  1 
their  disturbances  grow  monotonically .  Resolution  of  this  apparent  contradic 
tion  undoubtedly  lies  in  the  result  (30)  which  suggests  that  the  SVF  analysis 
by  restricting  attention  to  evolution  on  the  time  scale  t  =  0(0),  filters 
out  a  pulsating  mode.  In  fact,  Rogg  (1982)  has  reported  numerically 
determined  pulsations  for  flames  that  would  otherwise  be  candidates  for  SVF 
analysis. 

U.  Flames  Anchored  to  Burners 

We  now  turn  our  attention  to  gases,  which  necessarily  have  finite  Lewis 
numbers,  almost  invariably  lying  to  the  left  of  the  right  stability  boundary 
in  figure  5.  .  The  problem  is  to  find  a  mechanism  that  will  shift  this 

boundary  to  the  left ,  making  it  accessible  to  mixtures  of  practical  interest. 
Joulin  and  Clavin  have  shown  that  such  a  mechanism  is  the  distributed  heat 
loss  of  section  3.  ,  which  suggests  that  heat  loss  to  a  burner  anchoring 

the  flame  will  have  the  same  effect.  There  is  now  experimental  evidence 
■chat  burner  flames  can  indeed  pulsate. 

Consider  a  flame  anchored  to  the  so-called  porous-plug  hurner.  The 
mathematical  problem  to  be  solved  is 

3T/3t  +  3T/3x-V2T  =  -3Y/3t-3Y/3x  +  L_1V2  Y  =  DYe"9/T  with  V  =  DM"2,  (36) 

t  =  v  .  Y  —  L"13Y/3x  -  J  at  x  =  0,  T  bounded,  I  +  0  as  x  +  +«  (37) 

*s  s  • 

r.j  i3  the  prescribed  mass  flux  through  the  face  of  the  plug  at 
‘  r 

x  =  0;  whi'1®  T  and  J  are  the  prescribed  temperature  and  mas3-flux 
■  s  s 

fraction  there,  (ir.  practice,  cooling  coils  are  used  to  maintain  T 

s 

constant.)  The  physical  idea  underlying  the  boundary  condition  (37b) 
is  that  the  porous  surface  inhibits  the  flux  of  reaction  products,  so 
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that  the  flux  fractions  of  the  mixture  as  supplied  by  the  burner  are 

identical  to  the  mass  fractions  in  the  reservoir  supplying  the  burner. 

Analysis  of  the  steady  problem  in  the  limit  9  ■+  »  proceeds  as  in 

section  2.  for  the  unbounded  flame,  except  that  the  burnt-gas  temperature 

T.  is  the  fundamental  quantity  to  be  determined,  not  M  .  Integration  of  the 
o  r 

steady  version  of  equation  (36a)  from  x  *  0  to  ®  yields 

4 

T*  =  T  +  J  -  T  (38) 

s  s  s  D 


i.e.  the  heat  received  by  the  plug  in  terns  of  T^.  This  enables  us  to 
write  the  solution  in  the  form 

(j  -J1-^(T  +J  -T  )LeLx 

Js  S  S  S  D 


T  = 


f T  -J  +(T  +J  -T,  )eX 
1  b  s  s  s  b 


i 


Y  = 


for  x  >  x„,(39) 


where 

x*  =  £n[  J  / (T  +J  -T.  ) ] .  (UO) 

*  S  S  S  0 

Consistency  requires  x#  to  lie  between  0  and  “,  so  that  we  must  have 


T  <  T,  <  T  +  J  .  (Ul) 

s  b  s  s 

The  right  inequality  shows  that  Tj.  is  necessarily  positive,  i.e.  the 
plug  must  be  a  heat  sink.  Finally,  a  flame-sheet  analysis  gives 

M  =  /2lD  Tf  exn(-9/2T,  )/0J  ,  (U2) 

r  b  b  s 

from  which  T,  can  be  determined.  (The  result  (2.  )  is  recaptured  on 

D 

reolacing  J  with  Y„.)  Ferguson  and  Keck  have  made  satisfactory 

S  i. 

comparisons  between  experiment  and  a  theoretical  result  essentially 
equivalent  to  the  determination  (i*0)  of  the  stand-off  distance  x#  as 
a  function  of  the  injection  rate  M  •  (Note  that  is  also  used  in  making 
x,  nondimens ional . ) 
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The  inequalities  (hi)  define  limits  on  M  •  When  the  injection  rate  is 

r 

decreased  (increased)  beyond  its  limiting  value  a  surface  (remote)  flame  is 
obtained,  requiring  a  different  asymptotic  analysis.  We  shall  be  concerned 
only  vith  injection  rates  within  the  limits. 

5.  Stability  of  Burner  Flames 

In  considering  the  stability  of  the  solution  in  the  last  section  it  is^ 
natural  to  turn  to  a  HEP  analysis.  However,  such  an  analysis  requires  not 
only  L  to  be  sufficient  close  to  1  but  also  that  the  boundary  conditions 
permit  T  +  Y  to  be  constant  to  leading  order.  In  general,  the  conditions 
(37)  do  not  satisfy  the  second  requirement. 

One  way  out  of  the  dilemma  is  to  abandon  activation-energy  asymptotics 
and  adopt  a  suitable  modification  of  the  delta- function  model  used  in  the 
discussion  of  thermites.  (The  strength  of  the  delta  function  is  again 
proportional  to  exp(-6/2T#) . )  Such  an  approach  was  used  by  Margolis  who, 
by  means  of  a  numerical  investigation  of  the  complicated  dispersion 
relation  obtained  from  a  linear  stability  analysis,  was  the  first  to 
demonstrate  the  leftward  displacement  of  the  stability  boundary  alluded  to 
at  the  beginning  of  the  last  section.  He  also  carried  out  a  complete 
numerical  simulation  of  a  fuel-rich  hyrogen-oxygen  flame,  thereby  demon¬ 
strating  pulsations;  these  were  apparently  confirmed  by  experiments  performed 
at  Sandia-Livermore ,  although  there  is  no  published  account  of  them. 

The  NEF  requirement  that  T+Y  should  be  constant  to  leading  order  is  a 
sufficient  but  not  necessary  condition  for  the  flame  temperature  to  vary 
by  0(9 ~")  only.  Thus,  Buckmaster  (1982)  also  approached  the  question 
by  means  of  a  delta-function  model  but  then,  a  •posteriori,  justified  the 
model  through  activation-energy  asymptotics.  This  amounts  to  identifying 
the  circumstances  under  which  the  dispersion  relation  is  asymptotically 
meaningful,  and  for  which  there  is  a  flame  structure  linking  the  states 
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on  the  two  sides  of  the  flame  sheet  obtained  by  the  equivalent  jump  conditions. 
Some  care  is  necessary  because  of  the  sensitivity  of  the  solution  to  variations 
in  flame  temperature;  the  disturbance  field  is  0(l)  on  the  0-scale,  which 
entails  calculating  the  flame-temperature  perturbations  correct  to  0(9-^). 

In  turn,  this  entails  deriving  Jump  conditions  (from  the  flame-sheet  structure) 
correct  to  the  same  order.  Buckmaster  found 

S(  T)  =  -S(Y)  =  eTn/0,  (43) 

6(3T/3x)  =  -6(3Y/3x)  +  tJgW/0  =  -JgW+q/9  with  W  =  l+eT#1/2T^.  (44) 

Here  l  has  the  same  meaning  (4.  )  as  in  NEF  analysis,  e  is  the  small 

parameter  characterizing  the  size  of  the  disturbance,  the  flame  temperature 

is  T.  +  eT_,  /9 ,  and  q  is  a  quantity  (calculated  from  details  of  the  flame 

structure)  that  is  never  needed.  The  term  -  J  W  is  the  linearized  form 

s 

of  the  term  -Y  „W  appearing  in  the  Jump  condition  (20b),  with  J  replacing 
1  s 

V 

With  the  jump  conditions  in  hand,  it  is  a  straightforward  matter  to 
carry  out  the  analysis  of  unstable  disturbances  proportional  to  exp(at+iky). 
Circumstances  Justifying  the  delta-function  model  are  then  found  to  be 

9e"*X*  =  0(1)  (1*5) 

where  k  has  the  definition  (5*  )•  An  asymptotically  self-consistent 

dispersion  relation 

2<2(1-k2)+T(1+x)[  (l-02-lk2]  =  80x^e  KX*  (46) 

is  obtained,  the  corresponding  combustion  field  matching  the  structure 

used  to  deduce  the  Jump  conditions  (43),  (44);  here  l  and  0  have  the  definitions 

(U.  )  and  (21b)  with  in  place  of  Yf. 
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A  sinilar  dispersion  relation  (but  free  of  0  and  with  k  =  0)  was  derived 
by  Matkovsky  and  Olagunju  for  a  somewhat  artifical  burner  whose  boundary 
conditions  are  compatible  with  NEF  analysis.  They  do  not  discuss  the  full 
ramifications  of  their  results. 

What  emerges  from  Bucknaster 1 s  analysis  is  essentially  a  NEF.  To  get 
a  hint  of  this,  note  that  the  boundary  conditions  imply  that  the  disturbance 
satisfies 

Tx  +  Y1  =  L"13Y1/3x  at  x  =  0.  (47) 

Nov  the  disturbance  field  decays  rapidly  ahead  of  the  flame  sheet,  because 
Ax*  is  logarithmically  large  in  0.  It  follows  that  3Y^/3x  is  small  at  the 
plug,  small  enough  for  the  condition  (47)  and  the  near-equality  of  thermal 
ar.d  mass  diffusions  to  ensure  that  T^+Y  is,  at  most,  O(0-1)  throughout  the 
combustion  field. 

Finally,  it  should  not  be  thought  that  the  requirement  (45)  is  a 
constraint  on  x*.  Insofar  as  the  right  stability  boundary  is  concerned, 
the  relation  (46)  implies  that  as  x*  is  decreased,  Re(<)  increases 
(through  an  increase  in  the  frequency  of  pulsations)  so  as  to  keep  the 
term  on  the  right  balanced.  When  x*  is  0(l)  the  frequency  is  logarithmically 
large  in  0. 

The  displacement  of  the  right  stability  boundary  is  illustrated  by 
figure  3,  which  shows  how  the  point  k  «  0  on  it  varies  with  x*.  (Note  that 
l  -*•  2(l+/3)  as  x*  -*■  °°,  in  agreement  with  the  result  in  section  5.  for 
unbounded  flames.)  As  x*  Is  decreased  the  boundary  first  moves  to  the 
left,  by  an  amount  that  increases  with  0.  But  eventually  this  motion  is 
halted  and  the  boundary  moves  back  to  the  right. 
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6.  Pulsations  for  Rear-St agnation  Point  Flow. 

It  seems  likely  that  there  are  other  mechanisms  that  will  make  the 
pulsations  accessible.  One  that  has  been  suggested  by  theory  (but  not 
confirmed  experimentally)  is  negative  strain,  such  as  is  experienced  by 
a  flame  in  a  rear  stagnation-point  flow. 

For  moderate  Reynolds  numbers ,  flames  can  be  stabilized  behind  the 
closed  laminar  wake  at  the  rear  of  a  thin  plate  or  rod.  Figure  4  shows 
the  configuration,  which  Mikolaitis  &  Buckmaster  (1981)  treated  with  a 
NEF  formulation  based  on  the  equations 

(3/3t+ey3/3y-32/3y2)(T,h)  =  j£y3(0,T)/3y  (48) 

and  the  jump  conditions  (4.11--  •These  equations  are  the  unsteady  version 
of  equations  (4.^-0)  with  the  sign  of  €  changed  because  we  are  dealing  with 
a  rear  instead  of  a  front  stagnation  point.  The  problem  is  completed  by 
the  boundary  conditions 

T  =  T  ,  h  =  0  at  y  =  0,  (49) 

the  latter  corresponding  to  the  prescription  Y  =  Yf  at  y  =  0. 

The  steady  solution  can  be  written  in  a  closed  form  similar  to  that  in 
section  4.  for  a  front  stagnation  point.  Its  stability  to  one-dimensional 
disturbances  was  carefully  explored  using  a  combination  of  Galerkin's  method 
and  the  method  of  weighted  residuals.  The  results  are  illustrated  by  the  two 
curves  in  figure  5,  which  shows  variations  in  the  stand-off  distance  y#  of  the 
flame  with  (negative)  straining  rate.  All  responses  have  the  form  of  a 
backward  C,  so  that  there  is  a  maximum  straining  rate  beyond  which  the 


flame  must  blow  off . 
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_  2 

For  values  of  l  (=AY_/2T,  )  less  than  about  0.91  the  response  is  character- 

f  o 

ized  by  the  curve  (a) ;  the  lover  branch  is  stable  and  the  upper  branch  is 
unstable,  vith  a  real  eigenvalue  crossing  through  the  origin  as  the  turning 
point  T  is  traversed.  For  other  values  of  £,  as  characterized  by  the 
curve  (b) ,  part  of  the  lover  branch  is  also  unstable,  vith  a  complex 
conjugate  pair  of  eigenvalues  crossing  the  imaginary  axis  as  the  point  P 
(for  pulsations)  is  traversed.  This  raises  the  possibility  that,  for 
sufficiently  large  values  of  the  Levis  number,  blov-off  vill  in  practice 
be  preceded  by  pulsations.  (The  critical  value  of  l  is  quite  accessible.) 
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Figure  Captions 


7.1  Linear  stability  regions  for  thermite  flames,  with  admissible 
values  of  k  when  confined  to  insulated  circular  cylinder. 
Labels  to  the  right  (left)  of  the  neutral  stability  curve 
correspond  to  a  ■  4(2)  . 

7.2  Isotherms  of  spinning  thermite  flame,  with  +/-  denoting  the 
hot/cold  sides. 

7.3  Displacement  of  the  right  stability  boundary  in  figure  6.8 
due  to  the  presence  of  a  burner,  for  various  values  of  9  . 
The  purely  imaginary  numbers  are  values  of  a  . 

7.4  NEF  in  a  rear  stagnation-point  flow. 


7.5  Variation  of  stand-off  distance  y*  with  straining  rate  e 
for  NEFs  in  rear  stagnation-point  flows. 
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